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Abstract 

We have investigated theoretically the conductance of a Normal-Superconductor point-contact 
in the tunnel limit and analyzed the quantum interference effects originating from the scattering 
of quasiparticles by point-like defects. Analytical expressions for the oscillatory dependence of the 
conductance on the position of the defect are obtained for the defect situated either in the normal 
rS^ . metal, or in the superconductor. It is found that the amplitude of oscillations significantly increases 

when the applied bias approaches the gap energy of the superconductor. The spatial distribution 
of the order parameter near the surface in the presence of a defect is also obtained. 
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I. INTRODUCTION 



Electron scattering by single surface [1] and subsurface [2] defects results in an oscillatory 
dependence of the Scanning Tunnelling Microscope (STM) conductance G on the distance, 
r , between the contact and the defect. These oscillations originate from the interference of 
electron waves, which are scattered by the defect and reflected back by the contact. They 
have the same period (G ~ sin (2kpr + 5), kp is the Fermi wave vector) as the Friedel 
oscillations [3] of the local electron density of states in the vicinity of a scatterer. For 
subsurface point-like defects the oscillatory dependence of the conductance in a STM-like 
geometry has been investigated theoretically in Refs. [4-8]. 

Although defects below a metal surface can be 'visible' in STM data for up to ten in- 
teratomic distances [9, 10], the amplitude of the quantum oscillations in the conductance 
become very small with increasing defect depth. An effective way to enhance the STM 
sensitivity to such oscillation effects is to use a superconducting tip [11]. In Ref. [12] using a 
low-temperature STM with normal metal tungsten tips and superconducting niobium tips, 
the formation of electron standing waves near surface defects and step edges on a Au (111) 
surface have been observed. It was demonstrated that the amplitude of conductance oscil- 
lations is significantly enhanced when a superconducting tip is used, and when the applied 
bias \eV\ is close to the gap energy Ao of the superconductor. 

The investigation of various defects in superconductors with STM is of interest by it- 
self. For example, in Ref. [13] a bound state near a magnetic Mn adatom on the surface of 
superconducting Nb was observed by STM. The effect of single Zn defects on the supercon- 
ductivity in high-T c superconductors was investigated in Ref. [14], and the manifestation 
of d-wave symmetry of the order parameter was observed in the quasibound state near the 
defect. 

The listed reasons define the interest of theoretical investigations on the conductance of 
normal metal - superconductor (NS) tunnel contacts of small lateral size, in the vicinity 
of which a single defect is placed. The authors of Ref. [15] considered the conductance 
of a NS contact of finite size at low temperatures and for voltages \eV\ < A using the 
tunnelling Hamiltonian approximation. They found that, when the radius a of the contact 
is smaller then the Fermi wave length Xp, the conductance of a NS point-contact becomes 
G ns = (h/2e 2 ) G\ n ~ a 8 , where G nn is the conductance of the contact in the normal state 
[15]. This dependence is fundamentally different from the result of a quasiclassical theory 
[16], valid for a ^> \p. 

The conductivity of large (a ^> Xp) ballistic NS contacts in the presence of a 'planar 
defect' was investigated theoretically in several papers [21-24]. In these papers a planar 
NS structure and a 5-functional potential barrier, playing the role of the defect, have been 
considered, from which 'geometrical' resonances resulted due to combined Andreev and 
normal reflections. 

In order to describe the effect of isolated point-like defects in a superconductor on the STM 
conductance usually calculations of the local density of states n(r) are used (for a review, see 
[25]), where it is assumed that the conductance of the small tunnel contact is proportional 
to the local density of electron states. While for subsurface defects this assumption remains 
qualitatively valid, it does not permit a correct description of the details of the conductance 
oscillations because the bulk electron density of states around the defect is modified by 
reflection from the interface, r e S, and in the limit of zero tunnelling probability we have 
n(reS) = 0. In this case, the problem of electron transmission through the small NS 



tunnel junction in the presence of the defect should be considered. 

In this paper we present the results of a theoretical investigation of the conductance of 
a NS point contact (with a <C Xp ) in the tunnelling limit and we analyze the quantum 
interference effects originating from the scattering of quasiparticles by a point-like defect. 
Analytical expressions are obtained for the dependence of the conductance on the position 
of the defect and on the applied voltage, for the defect situated in the normal metal or in 
the superconductor. 



II. MODEL AND BASIC EQUATIONS 

Our model is presented in the Fig.l. The normal and superconducting half-spaces are 
separated by an infinitely thin dielectric interface, which has an orifice of radius a. The 
potential barrier in the plane of interface z = is taken to be a S— function, U (r) = 
Uof (p) S (z) , where p is the value of the radius vector p in the plane z — 0. The function 
/ (p) — > oo in all points of the plane except in the contact (p < a) , where / (p) — 1. In the 
point r a nonmagnetic defect described by a spherically symmetric potential D (|r — r |) 
is placed. A voltage V is applied between the two sides of the contact. We assume that 
the transmission probability \t\ of electrons through the barrier in the orifice is small (\t\ m 
t^kp/rrfUQ <C 1, m* is effective electron mass). In that case the applied voltage drops 
entirely over the barrier and the electric potential can be described by a step function, 
V (z) = V (—z) with V a constant. Based on the same reasoning we use a step function 
for the superconducting order parameter A (r) = A (r) (z). We consider the case of low 
temperatures and in the calculations take T = 0. At zero temperature a tunnel current flows 
through the contact for \eV\ > A . The applied bias is assumed to be small on the scale of 
the Debye frequency u D and the Fermi energy ep, \eV\ <C hu D <C e F . 

For definiteness we consider electron tunnelling from the normal half-space (z < 0) to 
the superconducting half-space (z > 0), i.e. eV > 0. In order to evaluate the total current 
through the contact, / (V), and the differential conductance, G (V) = dl (V) /dV, we should 
find the current density jk (r) of quasiparticles with momentum k at z > 0, formed by 
electrons transmitted through the contact. The current density jk (r) can be expressed in 
terms of the coefficients «k (r) and t>k (r) of the canonical Bogoliubov transformation [17, 18] 

j k (r) = 4 Im K(r)V<(r)/ F (E k ) - v k (r)Vv* k (r)f F (-E k )] , (1) 

where /p (E) is the Fermi function, which at T = is simply the unit step-function, /p (E) = 
(E) . The functions u k (r) and v k (r) satisfy to the Bogoliubov-de Gennes (BdG) equations 

[19] ' 



1,2 v 2 -e F + J D(|r-r |) 



2m 

h 2 



-V 2 -e F + D(\r-r \ 



Mk (r) + A (r) v k (r) = E k u k (r) , (2) 

v k (r) + A* (r) u k (r) = E k v k (r) . 



1m 
Eqs. (2) may be interpreted as wave equations for a two-component 'wave function', 
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V(z) = V 



z=0 



FIG. 1: Model of the contact. The point-like defect is situated in the normal half-space. The elec- 
tron trajectories in the normal metal and the trajectories of 'electron-like' and 'hole-like' excitations 
in the superconductor are shown schematically. 



of quasiparticles with energy E^. The conditions, which connect the vector ip k in the normal 
metal (V'nk) an d m the superconductor (ip s ^) &t the interface z = are 



^nk (P, 0) = ip s k (p, 0) = ip k (p, 0) 

2m* 



d - 9- 

T^sk [P, 0) - -Q^Pn-k (P, 0) 



ft 2 



^0/ (P) ^k (P, 0) 



(4) 
(5) 



The order parameter in the superconductor should be determined from the self-consistently 
condition 

A(r)= 7 Y, u k (v)vt(r)[l-2f F (E k )}, (6) 

A (z ^ +oo) ^ A , (7) 

where the constant Ao can be chosen real; 7 is the pair potential constant. It can be easily 
shown [17] that Eq. (1) combined with the self-consistently condition (6) automatically 
satisfies to the continuity equation 



divj^jk 



0. 



(8) 



The current-voltage characteristic / (V) of the contact in the presence of a defect can be 
found by means of integration of the current density jk (r) over the momentum k (within 
the energy interval Aq < E^ < eV ) and over a surface overlapping the contact in the 



superconducting half-space. For this surface we choose a half-sphere of large radius r ^> vq, £ 
(£ is the coherence length of the superconductor) centered at the contact r = 0. On this 
half-sphere we assume A (r) = Ao and hence E^ = \J £ k + Aq, where £ k = h 2 k 2 /2m* — ef 
is the kinetic energy measured from the Fermi level. The conductance G (V) of the contact 
(at T = 0) is given by 

oo 

G (V) = Anr e 2 N (0) J ^9 (z) J df k J ^0 (k z ) (rj k (r)) 6 (E k - eV) . (9) 

— oo 

where dfl and <if2 k are elements of solid angle in the real and momentum spaces, respectively, 
TV (0) is the density of states for one direction of spin. 



III. SOLUTION OF THE BOGOLIUBOV - DE GENNES EQUATION 

Generally, a self-consistent solution of Eqs. (2) can be found only numerically. Such 
solution must fulfil the condition of conservation of the total current / through any surface 
overlapping the contact, in spite of the spatial dependence of the order parameter. In order 
to simplify the task we will exploit the condition of a small barrier transparency and find 
an analytical solution of Eqs. (2) using the approximation of a constant order parameter 
A (r) = Ao© (z) . By means of this solution the coordinate dependence of A (r) can be found 
(see Appendix). 

In this section we generalize the method developed in the papers [4, 20]. We search the 
solutions of Eqs. (2) as an expansion into a series over the small transmission probability 
\t\ ~ 1/U , 

?k(r)=?ko(r) + 5ki00 + ---, (10) 

where ?/> k0 (r) satisfies the zero-boundary condition at z — 0, and -?/> kl (r) ~ 1/Uq- For the 

calculation of the current in leading approximation in the transmission coefficient (/ ~ 1/C/q ) 

it is enough to find the first correction ?/> kl (r). Substituting the expansion (10) into the 

boundary conditions (4), (5) we find that the function ip- kl (r) satisfies the condition of 

continuity at z — 0, and its value at z = +0 (in the superconducting half-space) is given by 

the relations 

h 2 d 

«ski(p,0) = ————-— M nko (p,0); Uski(p,0) = 0. (11) 

2m*U j (p) oz 

The boundary condition does not contain Andreev reflections, which appear in the next 
approximation in 1/Uq [30]. Thus, we will not consider Andreev resonances, which were 
analyzed in Refs. [21-24] for a one-dimensional model. 

The quasiparticle scattering by the defect will be taken into account by perturbation 
theory in the strength of the interaction with the defect. First, we find the solution of Eqs. 
(2) for the contact without defect. 

Let us consider an electron with energy i? k > A , which moves towards the interface 
from the normal metal. When D (r) = (the defect is absent) and 1/Uq = (the interface 
is impenetrable for electrons), in the normal half-space we have 

Mnko(r)=e^(e^-e-^), v nk0 (r) = 0, (12) 



where k = (x,k z ) , k z = fccos(i?), d is the angle between the vector k and the z axis, and x 
is the component of the wave vector parallel to the interface. 

Making use of the Fourier transform of the -*/> k (r) components over the coordinate p in 
the plane parallel to the interface, 



V>ki (p. *) 



dx / $ kl fx', z) e iVp , 



(13) 



and finding \I/ kl (x 7 , 0) from the simplified boundary condition (11), we find the solution of 
Eqs. (2) in the superconducting half-space 
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Mk^z 



f(p) 



fr(±) 



v/2 



m* 



/i 



ft 2 - 2 



x 



»n 



e F - 



1-Un 



2m* 
1 



ql/2 



:W£ 2 -A 2 
£ k 



2 V 1 + E k 



(14) 
(15) 

(16) 

(17) 
(18) 



£ (&z) — h 2 k z /im*Uo is the amplitude of electron wave after tunnelling through the ho- 
mogeneous barrier with a large Uq. The functions w k i (r) and t> k i (r) contain the sum 
of two solutions <p (r) of Eqs. (2), which correspond to 'electron-like' (k z > k z p = 
\^/lm* [ef — h 2 x 2 /2m*)) and 'hole-like' (k z ~' < k z p) quasiparticles having a positive z- 



component of the group velocity \ g = dE^/hdk. 

For a small radius of the contact (in the limit a — » 0) the function (16) takes the form [8] 



,(±) 



^ (r, k) 



(k^ a) 2 cos 9 



k±{E v 



\Jlvn* 

h 



e F ±\ E, 



h ( i\k {±) r), 

1/2 



A* 



(19) 
(20) 



Here, h\ (x) is the spherical Hankel function of the first kind. 

In the presence of the defect the functions Uki (r) and i>ki (r) can be found in first ap- 
proximation in the potential D (\r — tq\) of electron-impurity interaction by means of the 
Eqs. (2). 

1) If the defect is situated in the normal half-space the functions ii k i ( r ) an d fki ( r ) m the 
superconductor have the same form as Eqs. (14), (15) in which the amplitude t(k z ) must 
be replaced by the value 



t (k z ) = t (k z ) H -r 2 — gt (k) w nk0 (r ) h\ l> (kr ) 



,W 



(21) 



where g is the constant of the electron interaction with the defect 

g= fdrD(\r-r \). (22) 



In order to obtain Eq.(21) we assume that the characteristic radius of the scattering potential 
is much smaller than the Fermi wave length A^ (point defect). This condition permits taking 

the functions w n ko( r ) an d h\ (kr) outside the integral at the point r = r . The variations in 
the amplitudes of the 'wave functions' liki (r) and i>ki (r) result from the fact that the wave 
incident to the contact is a superposition of a plane wave and a spherical wave that comes 
from the scattering by the defect. 

2) If the defect is situated inside the superconductor, the additions A-u^i (r) and Ai^ (r) 
to the functions (14), (15) due to the defect scattering take the form 



Awn (r 



2nm*g 1 



dKe Mp- PQ ) (_L Mosiri (fc(+) z ) e ik ™ z [u u kl (r ) - *Wki (r )] 



— oo 

+—vo sin (k^z) e~ ik ' >z [u v kl (r ) - v u kl (r )] > ; (23) 

k z ) 

oo 

At; k i (r) = 2 ™* 9 y2 l _ u2 J dxe l ^P-P Q ) I-JL Vq sin (*.(+)*) e*™' [u Q u kl (r ) - ^o^ki (r )] 

— oo 

^u sin (k^z) (T ik< > )z [wo^ki (r ) - v a u kl (r )] \ . (24) 



fc.. 

It is known that the order parameter A (r) displays Friedel-like oscillations near a defect 
[26, 27] or a surface [28, 29]. The current through the tunnel contact / is defined by the 
average value of A (r), which coincides with Ao- In the Appendix we analyze the spatial 
dependence of A (r) near the surface of the superconductor, in the vicinity of which a non- 
magnetic defect is placed (at the distance less than the coherence length £ ). Figure 2 
illustrates the results of these calculations. An inhomogeneous spatial distribution of the 
order parameter is visible. We removed from the plot the region of radius A> (black circle) 
near the defect where Eq. (A9) is not valid. 

IV. CONDUCTANCE OF THE CONTACT 

By means of the solutions of the BdG equations, which have been obtained in previous 
section, we calculated the conductance G of the NS tunnel point contact. In linear approxi- 
mation in the electron-defect interaction constant g the conductance G can be presented as 
the sum of two terms, 

G (V, r„) = G 0ns (V) + AG osc (V, r ) , eV > A . (25) 

The first term, Go ns {V), in Eq. (25) is the conductance of the NS tunnel point contact in 
the absence of the defect 

Go n s (V) = Go nn — ^=j G 0nn = ., , , (26) 

\eVf - Al 97rnC/ o 




5 

xiXJln units) 



FIG. 2: Real space image of A (r) /Ao near the surface of the superconductor in the plane passing 
through the defect which has been obtained by using Eq. (A9), and the parameters zq = lOA^, 
£ = 10 4 X F , g = 4ir. 

where Go nn is the conductance of a contact between normal metals, which is multiplied by 
the normalized density of states of the superconductor at E — eV in Eq. (26). The second 
term describes the oscillatory dependence of the conductance on the distance between the 
contact and the defect. 

If the defect is situated in the normal metal half-space AG osc (V, r ) is given by 



AG osc (V, r ) = -G 0ns (V) — g — (k F z ) j 1 (k F r ) y 1 (k F r ) , 

it \r 



(27) 



where 



2nm*k I 



(28) 



is the dimensionless electron-defect interaction constant, ji(x) and yi(x) are the spherical 
Bessel functions of the first and the second kind [31], and X F = h/ \/2m*e F . In Fig. 3 
dependencies of AG osc (V, tq) on the distance p are shown for two values of the bias eV, one 
of which is very close to the gap energy (eV/Ao = 1.1), and the second one is eV = 2Ao- The 
figure illustrates the increasing amplitude of the conductance oscillations near eV ~ Ao- 

For the defect in the superconducting half-space the oscillatory part of the conductance 
consists of two terms 



12_ / X \ 

AG osc (V, r ) = -G 0ns (V) —g — (k F z ) 2 V] ^ a (eV) ji (k a r ) y t (k a r ) 



(29) 



where 



1>i 



v 



V2 



m 



h 



e F ± 



^f(^Vf 



nl/2 



A 2 



(30) 
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FIG. 3: Dependence of the normalized oscillatory part of the conductance AG OSC /Go ns , Eq. (27), 
on the distance p between the defect and the contact axis for two values of the applied voltage. 
The defect is situated in the normal metal at a depth zq = 5Xp. The dimensionless constant of 
interaction is taken as g = 0.01. 

In Eqs. (26)-(29) we neglected all small terms of the order of A /ep and eV/ep . Nevertheless 
we kept the second term in square brackets in the formula for k± (see, Eq.(30)) because for 
a relatively large r , (\/(eV) 2 — A^/sp^ro/Air) ~ 1, the phase shift of the oscillations may 
be important. In Fig. 4 we show the difference between the dependencies of the normalized 
oscillatory parts of the conductance AG OSC /Go ns on the distance p for a contact between 
normal metals (A = 0) and for a NS contact. An observable shift of the conductance 
oscillations results from the voltage dependence of wave vectors k± (30). 

V. CONCLUSION 



Thus, we have analyzed the conductance G of a tunnel NS point contact with a radius 
a smaller than the Fermi wave length Xp, at low temperatures (T = 0) and for applied bias 
eV larger than the gap energy of the superconductor Ao- The effect of quantum interference 
of quasiparticles scattered by a single defect situated in the vicinity of the contact has 
been taken into account. We have shown that in leading approximation in the parameters 
eV/ep <C 1, A /ep <C 1 the conductance of a small NS contact is G 0ns = G 0nn N s (eV), 
Eq. (26), i.e., the product of the conductance of the same contact between normal metals, 
Gonn ~ a 4 , and the normalized density of states of the superconductor N s (eV), similar as 
for a planar tunnel contact. Although such result is not unexpected and has been confirmed 
by experiment [11], for a contact of radius a < Xp it was not obvious and it is first obtained 
in this paper. 

If the defect is situated in the normal metal the conductance displays oscillations, the 
period of with is defined by the Fermi wave vector, AG osc (V, r ) ~ sin 2kpr at kpr ^> 1 
(Eq. (27), Fig. 3), as for a contact between normal metals [4]. In this case the defect plays the 
role of an additional 'barrier' between the normal and superconducting metals and results 
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FIG. 4: The dependence of the oscillatory parts of the conductance AG osc /Gq (29) on the distance 
p between the defect and contact axis for the contact between normal metals (AG^lc /Gonn) an d 
a NS contact (AGdsc /GWs). The defect is situated in the right metal (the superconductor) at a 
depth WX F ; eV/A = 5; g = 0.01. 

in oscillations of the transmission coefficient. The underlying principle here is similar to 
resonance transmission through a two-barrier system. 

In the superconductor the electron wave incident on the contact from the normal metal is 
transformed into a superposition of 'electron-like' and 'hole-like' quasiparticles. In the case 
of location of the defect in the superconducting half-space quantum interference takes place 
between partial waves transmitted and those scattered by the defect, for both types of quasi- 
particles independently (Eq. (29)). Although the difference between wave vectors k^ (eV) 
of 'electrons' and 'holes' is small the shift (k^ — fc( - )) r between the two oscillations should 
be observable (Fig. 4). 



Appendix: Oscillations of the order parameter near the surface in the presence of 
a defect. 



When calculating the conductance to first order in the transmission probability we should 
know the order parameter A (r) in the limit of a nontransparent interface (surface), Uq — *■ oo. 
According to Ref. 



A*(r) 



7 T 



OC 



^ F+(r,v)Q(u> D -u>) 



(Al) 



where u = ttT (2n + 1) are the Matsubara frequencies. The Fourier components G w (r, r') 
and F+ (r, r) of Green's functions satisfy the Gor'kov equations, which in the absence of a 



10 



defect potential have the form 



iuj 



h 2 v 2 

2m* 

h 2 v 2 



e F )G U} h,r') + A(r)F" 



(iu + ^ + £*■) F+ (r, r') + A* (r) G u (r, r') 



6 (r - r') 



0. 



(A2) 



For a homogeneous superconductor A (r) = An = const, and the solutions G w (r, r') 



Gi 0) (r - r') and F+ (r, r') 



-(o) 



r — r') of Eqs.(31) can be found to be 



a (0) 



p+(o) 



7iN(0) 

COS KpT 

KpT 

ttN (0) A* sin k F r 



IUJ 



'A 2 + u 2 



sin kpr 



exp 



Vph 



y/*Z + «> 2 ) ,(A3) 



Aj 



LO 



2 k F r 



exp 



Vph 



A 2 



w 



(A4) 



where r = |r — r'| , vf is the Fermi velocity, uo <C ep. For the semi-infinite superconducting 
half-space any component of the matrix Green function 

F^\. 
/■r'Mr.r') _(V'''- ' (A5) 



Gi s) (r, r' 



C 






^i(r',r) 



can be written as 

G« (r, r') = Gi°) (r - r') - G^ (r - r>) , (A6) 

where r' = (x 1 , y', —z') . Equation (A6) is exact and it provides the zero value of A (r) at the 
surface z — 0. The fact that the order parameter vanishes at the nontransparent interface 
can by seen from Eq.(6). 

The Green's function for the superconducting half-space in the presence of the point 
defect can be found from the Dyson equation 



G,,(r,r f ) = GW 



r, r 



J dv"G^ (r, r") D (|r"-r |) r s G w (r", r') 



(A7) 



where t% is the Pauli matrix. Making use of the small radius of the defect potential in the 
first order approximation in the interaction constant g (22) we obtain 

(A8) 





F+(r,r) 


= F^ (r, r' 


) + 






^ +(s) (r 


r )C#>(ro,] 


0+^i(r , 


r) F^ 


(r , 






As a first step for the self-consistent solution, the functions Gi, (r — r') (31) and 

Fa, (r — r') (31) may be used. At T — *■ the summation over Matsubara frequencies 
in Eq.(Al) can be replaced by an integration. Substituting the Eqs.(31), (31) into Eq. (A6) 
and using Eq.(A8) we find the space distribution of the order parameter (Al) in the next 
(after A = A = const.) approximation. 



1 



4^ ln " 



A(r) = 
2ui D 



sin2k F z _ x f2cu D 
A < 1 — In 



2k pz 



An 



A r 



sm2k F s f 2irs uop 



sin k F (so + s o' 



7C 



o/ 2(kps 

K (S + go) C^d 



sin2A; F so „ (2tcs uo d 



(A9) 



11 



Here 

arshb 



K (a; b)= dte- ach \ (A10) 

o 
s = |r — r | ; s = |r — r | , and £ = hvp/n^Q is the coherence length. At ab ^> 1, 
/C (a; b) ~ A" ( a )> the modified Bessel function [31]. The Eq.(A9) is valid at distances from 
the defect larger than the characteristic radius of the potential D (|r — r |) . The correction 
to the constant value of the order parameter A decreases at small distances r <C £ from the 
surface or the defect according to a power law, and vanishes exponentially (~ e~ 27rr ^°) at 
larger distances r ^> £ . A grey-scale plot of A (r) obtained by means of Eq.(A9) is presented 
in Fig. 2. In the plot we used an unrealistically large value of the constant g in order to 
show the influence on the order parameter the defect and the surface in the same plot. For 
realistic values g ~ 0.01 the spatial oscillations of A (r) resulting from the scattering by the 
defect have a much smaller amplitude than the second term in the braces of Eq.(A9). The 
matching procedure can be continued when we put A (r) of Eq.(A9) into Gor'kov's equations 
(Eqs.(31)) or BdG equations (2). Unfortunately, starting with this step the solutions may 
be obtained only numerically. 
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